Universal finite-size scaling analysis of Ising models with long-range interactions at 
the upper critical dimensionality: Isotropic case 
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We investigate a two-dimensional Ising model with long-range interactions that emerge from a 
generalization of the magnetic dipolar interaction in spin systems with in-plane spin orientation. 
This interaction is, in general, anisotropic whereby in the present work we focus on the isotropic 
case for which the model is found to be at its upper critical dimensionality. To investigate the critical 
behavior the temperature and field dependence of several quantities are studied by means of Monte 
Carlo simulations. On the basis of the Privman-Fisher hypothesis and results of the renormalization 
group the numerical data are analyzed in the framework of a finite-size scaling analysis and compared 
to finite-size scaling functions derived from a Ginzburg-Landau- Wilson model in zero mode (mean- 
field) approximation. The obtained excellent agreement suggests that at least in the present case the 
concept of universal finite-size scaling functions can be extended to the upper critical dimensionality. 
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I. INTRODUCTION 

In the last decade, spin models with long-range inter- 
actions were the subject of several extensive Monte Carlo 
studies. Utilizing an efficient cluster algorithm pj these 
studies were addressed to the verification of some un- 
proved predictions on the critical behavior of spin mod- 
els with algebraically decaying long-range interactions 
0- Furthermore, the crossover from Ising- like to clas- 
sical critical behavior was investigated [3, |3| and first 
numerical results on the critical behavior of the dipolar 
in-plane Ising (DIPI) model were obtained @ . This two- 
dimensional model displays a strongly anisotropic phase 
transition, i. e. , the correlation lengths in direction par- 
allel and perpendicular to spin orientation diverge in the 
infinite system (let t > 0) as [5| 
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at the critical point, where both £|| ^ £j_ and v\\ ^ v±, 
and t = (T — T c )/T c denotes the reduced temperature. 
Except the anisotropy exponent 9 — v\\/v±_ neither any 
numerical estimates of the critical exponents exist for the 
DIPI model, nor is it clear whether the model exhibits 
Lifshitz type critical behavior [|| as it is observed, e. g. , 
in the anisotropic next nearest neighbor Ising (ANNNI) 
model 



To address such questions in a broader context we 
present a two-dimensional long-range spin model that 
arises from generalizing the DIPI model by introducing 
an additional parameter. Assuming Lii x L± geometry 
and periodic boundary conditions this model is described 



Figure 1: The model (Eqs. @) exhibits a ferromagnetic 
ground state in the region denoted by "fm". The dashed lines 
u\i — — 2u;± > and — 2ixJ| = u>± > 0, respectively, and 
the origin wii = lu± = also represent ferromagnetic ground 
states. All other regions correspond to modulated spin con- 
figurations (see text). 



by the Hamiltonian 



(2a) 



with magnetic Ising spin variables Sj = ±1, the spin-spin 
distance vector r^ , and an external field B. The pair 
coupling J(r) is given by 



J(r) = JS\ r \ A + 



(2b) 
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and consists of both a ferromagnetic short-range nearest 
neighbor exchange coupling with the coupling constant 
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J > and a long-range contribution, where ru and r± 
are the components of the vector r parallel and perpen- 
dicular to spin orientation. Using this general form of the 
pair coupling J(r), several well known spin systems can 
be mapped onto this model by making an explicit choice 
of the parameters ton, tu±, and J. With u>u = —2lj± > 0, 
and for symmetry reasons — 2u>\\ = lo±_ > 0, and J > the 
DIPI model is recovered, and if u>n = oj± > and J = 
Eq. I|2a[) corresponds to an Ising model with an isotropic 
ferromagnetic long-range interaction algebraically decay- 
ing as J(r) oc | r | 3 . Another special case is the dipolar 
Ising model with perpendicular spin orientation jlOl [llj 
that can be obtained for uu = ui± < and J > 0. 

Figure shows the ground state phase diagram of the 
model whereby we took into account four different ground 
state spin configurations: the ferromagnetic state (fm) 
where all spins point to the same direction, the totally 
antiferromagnetic state that is referred to as checker- 
board state (check), and commensurate stripe domain 
states with a domain wall orientation parallel (mod ||) 
and perpendicular (mod _L) to spin orientation and the 
periods -/Vn j_. Dependent on the values of the quotients 
ui\\/J and uj±/J all considered spin configurations were 
found as stable ground states. Due to symmetry reasons 
the arrangement of the corresponding phases in Fig. ^ is 
symmetric with respect to the line Wii = u)j_. The periods 
-/Vim of the stripe domain states diverge when approach- 
ing the dashed lines in Fig. ^ 

The region in parameter space where the model dis- 
plays a ferromagnetic ground state is of particular inter- 
est to us. In that region the observed phase transitions 
are isotropic (anisotropic) when u>u = loj_ (cjii =/= u>±) 
whereby in the present work we draw our attention to 
the isotropic ferromagnetic long-range case 

W|| = cj_l > 0, J = (3) 

before we turn to the case of an anisotropic pair cou- 
pling |l2^. It is known that the upper critical dimen- 
sion of long-range spin models with ferromagnetic inter- 
actions decaying as J(r) oc |r|~( d+cr ) is given by d u = 2a 
|13| . Comparison with the pair coupling J(r) defined in 
Eq. i|2b|) yields in the isotropic case, Eq. ©,(7 = 1 for a 
two-dimensional system and consequently d = d u = 2. 

So in this paper we investigate the critical behavior of 
the model at its borderline dimensionality d u by means 
of Monte Carlo (MC) simulations and finite-size scaling 
methods. For that purpose in the following Section [Q] 
the finite-size scaling form of the free energy density is 
discussed and the finite-size scaling relations of the con- 
sidered quantities are derived as they are used for the 
finite-size scaling analysis. These relations define finite- 
size scaling functions for which in Section lrTTl wc evaluate 
analytical expressions in the framework of the so-called 
zero mode theory that is based on the Ginzburg-Landau- 
Wilson (GLW) model. In the last Section II VI the zero 
mode results are compared to numerical data within a 
finite-size scaling analysis. 



II. FINITE-SIZE SCALING RELATIONS 

To study the critical properties of the model in the 
isotropic long-range case we have carried out a finite- 
size scaling analysis of MC data. This analysis requires 
the finite-size scaling relations of the quantities that were 
considered in the simulations. 

Via a renormalization group approach, Luijten and 
BlSte derived the scaling form of the free energy den- 
sity of 0(n) spin models with ferromagnetic long-range 
interactions decaying as J(r) oc |r| - ( rf+cr ). At the upper 
critical dimension, that is given by d u — 2 for a = 1 
(see SecQJl, the singular part of the reduced free energy 
density was found to scale as (n = 1) 

f s (t, h; L) ~ L- 2 /(u-*LIn*(^) [t - v' L' 1 , 
Lhni(^r)h) (4) 

with the reduced temperature t, the reduced external 
field h = (3B where (3 = l/(fcBT) denotes the inverse tem- 
perature (we set fce = l throughout this paper), and the 
so-called dangerous irrelevant variable u [T-U fl5| . Note 
that we rewrote the formula given in Ref. |2j in terms of 
the parameters v' and L' , where L' can be regarded as 
a reference length that fixes the length scale in the log- 
arithms (see also Ref. ^^), and we absorbed constant 
factors into /. The symbol ~ means "asymptotically 
equal" and, unless stated otherwise, refers to the limit 
(i, h, L) (0, 0, oo) with tL ln^ (L) and hhi In' (L) fixed 
(cf. Eq. 0). 

Proceeding from Eq. we adapt the Privman-Fisher 
hypothesis |Tj| and propose the finite-size scaling form of 
the singular part of the reduced free energy density 

f s (t,h;L)~L- 2 Y(x TS ,y IS ) (5) 

with the universal finite-size scaling (UFSS) function 
Y(x,y). The arguments of this function correspond to 
the temperature scaling variable 

x ts = C 1 iL]n*{±) (6a) 

with the shifted reduced temperature 

t = t-vL- 1 \n-i(±), (6b) 

and the field scaling variable 

y Ig = C 2 hLhn*(%), (6c) 

whereby C\ and C2 are nonuniversal metric factors. Let 
us note that we have replaced the constants v' and L' a 
which arise from the renormalization group by the con- 
stants v and L that will be used as fit parameters in the 
finite-size scaling analysis. 

It is also important to point out that the terms that 
result from v and Lq in the temperature scaling variable 
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x Ig and the field scaling variable y vg are merely correc- 
tions since they do not contribute to the leading orders 
in the expansions 

tLlJ(^) L =°° tL\J{L) +vO (In-? (L)) 

+tL\n(L )O(ln~*(L)) (7a) 



+hL* ]n(Lo)0(kr*(L)). (7b) 



However, due to the slow convergence of the logarithms 
appearing at the upper critical dimension these correc- 
tions are substantial for the quality of the data collapse 
in the finite-size scaling analysis, as will be discussed in 
Sec. H3 (see also Sec. IV(B) in ||, and [3). 

In the following we derive the finite-size scaling forms 
of the quantities that were considered in the simulations. 
Let 



(8) 



denote the average of the spin variables Sj , these quanti- 
ties are the magnetization m(t, h; L) = (s) and the sus- 
ceptibility x(t, h; L) = (3L d ((s 2 ) — (s) 2 ), whose finite-size 
scaling forms can be obtained by taking derivatives of the 
singular part of the reduced free energy density, Eq. JSJ, 
according to 

d 

m(t,h;L) = -—f s (t,h;L) 

~ C 2 L-hni(±)Y m (x Ig ,y Ig ), (9a) 

p- l X {t,h;L) = -^Lf,(t,h;L) 

~ ClL\J(±)Y x (x rg ,y Ig ). (9b) 

We also consider the dimensionless Binder cumulant 0] 
U(t;L) — 1 — (s 4 )/(3(s 2 ) 2 ) that can be evaluated from the 
susceptibility \{t, h; L) and the nonlinear susceptibility 
X {nl) (t,h;L) using the identity [iH ] 



U(t;L) 



X {nl) (t,0;L) 



(10) 



3/3LVM;£)' 
where x^(*> ^ L) is given by 

(3~ 3 X (al) (t,h;L) = -J^f,(t,h;L) 

~ C*L 4 ln(±)Y x w(x Ig ,y Ig ).(ll) 

Hence, the Binder cumulant scales as 

U(t;L) ~ Yu(x Ig ), (12) 

where Y v (x) = -F x („d (x, 0)/(3F x 2 (a;, 0)). Since the en- 
semble averages (s 2m+1 ) with m € No vanish in the ab- 
sence of an external field, we, in addition, analyze the 



magnetization m a b s (i, h; L) = (\s\) and the susceptibil- 
ity X a bs(t,h;L) = (3L d ((s 2 ) - (|s|) 2 ). It is understood 
that these quantities also fulfill the finite-size scaling 
forms of m(t,h;L) and x{t,h\L), Eqs. (J5J), respectively 
with the corresponding UFSS functions Y m ^ bs (x,y) and 



III. MEAN-FIELD THEORY 

An appropriate description of spin systems of dimen- 
sions LxLx---xL = L d with mean-field-like (classical) 
critical behavior can be achieved by the mean-field the- 
ory employ ed b y Brezin and Zinn- Justin |2£j and Rud- 
nick et al. (2lj | . This theory, also known as zero mode 
approximation, yields in contrast to conventional mean- 
field theories a rounded transition for finite systems. In 
the thermodynamic limit L — > oo the usual power laws 
with the expected mean-field values of the critical expo- 
nents can be recovered. 

In the following, this theory is reviewed and used to 
evaluate analytical expressions for the finite-size scaling 
functions defined in the preceding section in order to com- 
pare them to the numerical data, as it is demonstrated 
in Sec. IIVI The basis of this evaluation is the reduced 
GLW Hamiltonian in momentum space that corresponds 
to the underlying spin system with a long-range interac- 
tion J(r) cx \r\~ i - d+,T \ It is given by (see, e.g. , Ref. 0) 



H 



-41 ^2^2^2fk 1 (pk 2 fk 3 'p-k 1 -k 2 -k 3 ] 

ki k 2 k 3 / 



(13) 



with the temperature-like parameter r oc T — T c that 
measures the deviation of the temperature from the 
mean- field critical temperature T c m , the reduced exter- 
nal field h, and the dangerous irrelevant variable u > 0. 
Each sum in Eq. 1)13(1 runs for each component kj of k 
over integer multiples of 2ir/L up to a momentum space 
cutoff k\ = 7r / 'a (\kj \ < k\) with the lattice constant a. 

The essential step of the zero mode approximation is 
the neglection of all modes except the zero mode tp = 
ipo in Eq. l(T3]l. This leads to the reduced zero mode 
Hamiltonian (see, e.g. , Ref. |2lj ') 



(14) 



W fo>) - L*{-tp' + - hip 
with the corresponding partition function 

2„ = / dLpe- n °^\ 



(15) 



So the normalized probability distribution of the zero 
mode is given by 



1 

— ( 

Z 



-Wo(<p) 



(16) 
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and it can be used to evaluate averages of the form 

/>oo 

(g(<p)) = / dtpg(<p)V (<p). (17) 



A further central quantity is the reduced free energy den- 
sity that is given by f n = — L~ d h\(Zo) within the zero 
mode approximation |l9l |20| . Using this expression and 
the zero mode partition function defined in Eq. i|15|) , the 
rescaling (see, e. g. , Ref. poT ] 1 ) 



ip 



(uL d ) 



(18) 



immediately yields the zero mode finite-size scaling form 
of Jo- It reads 



f (r,h;L) = L d f {x mf ,y m f) + c(L) 



(19) 



with the mean-field temperature and field scaling vari- 
ables 



x m f = ru 2 L'- 



Vmf 



1 — I J- — 
hu 'L'. 



(20) 



and an additive term c(L) that is without significance in 
the following since it is absent after taking derivatives of 
/o(r, h; L) with respect to r or h. Instead we focus on the 
finite-size scaling function 



fo{x,y) = - In 



dip e 



-(§ v 2 +^v 4 -av) 



(21) 



from which, as seen in the preceding section, the finite- 
size scaling functions of other quantities like the magne- 
tization m and the susceptibility x follow. 
The asymptotics of this function are given by 



/o(*,0) 

/oM) 

/o(0,tf) 



1, (\x\ 



81 



32,' 12/13 



(22a) 
(22b) 
(22c) 



Following the convention suggested in Ref. |23J, the nor- 
malized finite-size scaling function Y mi (x,y) of the re- 
duced free energy density should be defined such that 
their asymptotics read 



Y mi {x,0) 
Y mt (0,y) 



X—> — OQ 



\y\ 



(23a) 
(23b) 



instead of the leading orders in Eqs. I|22l) . This require- 
ment fixes some arbitrariness of the reduced free energy 
density finite-size scaling function and can be fulfilled 
by a rescaling of the parameters r and u, which are the 
only phenomenological quantities entering the reduced 
zero mode Hamiltonian, Eq. 114(1 . Replacing x and y in 
Eqs. (|22|l explicitly with the mean-field scaling variables 
x m { and y m { from Eqs. I|20|l . the rescaling 



3\/3 



-r, 



81 
32 ? 



leads to a cancellation of the corresponding coefficients 
of — xtrf and — \y m f | a and one obtains the desired asymp- 
totics. It is important to note that we do not alter the 
definitions of x m f and y m { due to this rescaling, but the 
scaling function itself. Starting from the zero mode par- 
tition function, Eq. (|15J) . with the rescaled parameters r 
and u, we, after the procedure discussed above, finally 
end up with 



F mf (x, y) = - ln(S (x, y) + E (x, -y)) (25) 



where 



E m (ar,s/) 



d<p ip r ' 



(26) 



instead of fo(x,y) given in Eq. I|21|) . 

Having defined the reduced free energy density finite- 
size scaling function, in the following the finite-size scal- 
ing forms and the corresponding zero mode finite-size 
scaling functions of the quantities considered in the pre- 
ceding section will be derived. Since the quantities 
w a b s (i, h; L) and Xabs(t, h; L) cannot be evaluated by tak- 
ing derivatives of the reduced free energy density, we in- 
stead make use of the average defined in Eq. (|17f) . Set- 
ting here g(ip) = ip m and g(ip) — \ip\ rn , respectively, the 
rescaling of the parameters r and u [Eq. 124|) ] and the 
rescaling of the zero mode (Eq. (|18fl ) immediately yields 
the finite-size scaling forms of these averages. They read 



( 1 iLt^ (m) ' m '(^,!/ m f), (27a) 
(^ d )~^n ( b" l) ' mf (^mf,y mf ) (27b) 



with the finite-size scaling functions 

Y^ f (x,y) = + (-lrSnjx-y) 



E (x,y) + E (x,-y) 
E m (x,y) + E m (x,-y) 
S (x,y) + E Q (x, -y) ' 



(28b) 



Since, as discussed in Ref. [2(j, the zero mode ip is related 
to the order parameter field </>(r) in real space via 



(p ■ 



1 L 

i=i 



(29) 



and consequently corresponds to the average order pa- 
rameter per volume, one immediately obtains the finite- 
size scaling forms of the quantities defined in Sec. [H] 
using Eqs. f2*7| . Due to this correspondence, the mag- 
netization m(r, h; L) and the susceptibility x(r, h; L) are 
given by 

m(r,h;L) = (tp)o (30a) 
p^xfahiL) = L d ((ip 2 ) - (ip)*) (30b) 



(24) 

and therefore, according to Eq. (|27ap . scale as 
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with the two-dimensional lattice sum [2j 



m(r,h;L) = u^L^Y™* (x mf , y ml ) (31a) 
P^X^hiL) = u-?LiYf{x mt ,y m i) (31b) 

with the zero mode finite-size scaling functions Y™ (x, y) 
and Y™ l (x,y). The finite-size scaling forms of the 
quantities 



m ahs (r,h;L) 
/3" 1 Xabs(r,/i;L) 



(M)o) 



(32a) 
(32b) 



are identical to Eqs. i|31|) . respectively, with the cor- 
responding finite-size scaling functions Y^f {x,y) and 
Y^ e (x,y). A further quantity of interest is the dimen- 
sionless Binder cumulant 



U(r;L)= 1- 



<^ 4 



3(^ 



2\g 



h=0 



for which within the zero mode approximation one ob- 
tains the scaling form 



U(r;L) 



(X m f) 



(34) 



with the finite-size scaling function Y f mf (x). 

All zero mode finite-size scaling functions Y™ (x,y) 
can be expressed as combinations of the functions 
E m (x,y) defined in Eq. JSBJ with m = {0, 1,2,4}. Since 
to analyze the critical behavior, as it is done in the MC 
simulations, either the temperature scaling variable or 
the field scaling variable is kept at its critical point value, 
it is sufficient to evaluate S m (x,0) and S m (0, y) for the 
pertinent values of m. Analytical expressions for the 
needed functions can be found in the Appendix. 



IV. MONTE CARLO RESULTS 

In this section we present the results of MC simula- 
tions that were carried out for the model [Eq. J5J)] with 
(u)\\,w±,J) = (1,1,0). We have used the Wolff cluster 
algorithm p3 | for long-range spin models proposed by 
Luijten and Blote 0. To study the properties of the 
model in the presence of an external magnetic field, a 
histogram reweighting technique [25j was used. 

In the simulations quadratic spin systems with L := 
L|| = L± = {32,64,128,256,512,1024} were considered 
and we started with recording the magnetization rn a bs, 
the susceptibilities x and Xabs, and the Binder cumu- 
lant U at zero field for various temperatures close to the 
corresponding mean- field critical temperature T c mf . This 
temperature is given by T c mf = J(0) (see, e. g. , Ref. |2(|) 
where J(k) denotes the Fourier transform of the pair cou- 
pling Eq. (|2b|> . Setting the coupling constant to J = the 
evaluation of this expression yields, the divergent term at 
r = is excluded, 



e(s) = A 1 -x(s) (C(*,|)-CM)) 



(36) 



where £(s, a) denotes the generalized Riemann Zeta func- 
tion. Setting cjii = lo±_ = 1, the mean-field critical tem- 
perature takes the value 



9.0336. 



(37) 



According to the finite-size scaling relations listed in 
Sec. [H] the UFSS functions were evaluated from the MC 
data and plotted against the temperature scaling vari- 
able x rg (see Fig. EJ. The data collapse was achieved by 
adjusting the critical temperature T c = T c (oo) and the 
parameters v and Lq in the following way. First we de- 
termined L from the requirement that the maximum of 

(33) 

the scaled susceptibility Txabs(i, 0; L)L 1 In 2 (-£-) col- 
lapses for different L [Fig. |2Jb] , as this peak height is in- 
dependent of T c and v. After that we adjusted T c and v 
until the scaled cumulant U(t;L) vs. x rg [Eqs. (|12l6a|) ] 
collapses [Fig. |5JI] and fits the well known critical value 
ill 



Yu(0) = 1 



r(|) 4 

24tt 2 



= 0.27052. . 



(38) 



at x rg =0. Finally, in order to compare the numerical 
data to the zero mode finite-size scaling functions listed 
in Sec. IIII1 these functions were fitted to the numerical 
data by tuning the nonuniversal metric factors C\ and 
C*2 • The values of all parameters as they were determined 
from this analysis are listed in Table Q] The whole data 
analysis was done using Fsscale |28j . 

In addition to the temperature dependence we also 
studied the dependence of the quantities m, m a b s , X, and 
X a b s on an external field at the temperature 



T c (L) = T c (oo)(l + V L- 1 ln-3(^)) 



(39) 



that can be regarded as an effective critical temperature 
of the finite system (see, e.g. , Ref. ^^)- Note that 
t c {L) — T c (L)/T c (oo) — 1 corresponds to the value of t 
for which the shifted reduced temperature t as defined in 
Eq. (|6b|) and consequently the temperature scaling vari- 
able x rg [Eq. iJBaJ] vanishes, as i = t — t c (L). Due to this 
choice the numerical data obtained at this temperature 
and nonzero fields might then be compared to the corre- 
sponding zero mode finite-size scaling functions 5^ mt (0, y) 
within a finite-size scaling plot. Therefore we stored for 
each of the system lengths L — {16,32,64,128,256} a 



Tc(oo) 


L 


V 


Ci 


c 2 


8.0302(3) 


3.0(2) 


1.16(2) 


0.735(10) 


0.92(1) 



E 



1, 

= 2^11 



+ ua)6(§) (35) 



Tabie I: Values of the fit parameters used throughout this 
work. 
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^3 



■ L = 32 

* L = 64 
L= 128 

■ L = 256 

* L = 512 

* L = 1024 

- Y mi (x ,0) 

v rg' ' 




0.5 r 




■ L = 32 

* L = 64 
L= 128 

■ L = 256 

* L = 512 

* L = 1024 
- Y M (x ,0) 



C 1 [fLln 1/6 (L/V- vlnl/2 ( L/L o)] 



C, [tL\n\LIL)-v\n m {LIL)-\ 



-A 



-J 
© 




0.1 L 







■ L = 32 

* L = 64 
L= 128 

■ L = 256 

* L = 512 

* L = 1024 

— y mf (x ,o) 

X V rg' ' 



■ L = 32 
+ L = 64 

L= 128 

■ L = 256 
' L = 512 

♦ L = 1024 
- Yf{x ) 

U v rg y 



1/6 



1/2, 



x =CJtL\D. (L/L Q )-v]n (UL Q )] 




x ^[fLln (L/L Q )- vln" (L/L Q )] 



Figure 2: Finite-size scaling plot of the magnetization m a bs(£,0;L) (a), the susceptibility Xabs(i,0;L) (b), the susceptibility 
x(£, 0; L) (c), and the Binder cumulant U (t; L) (d), versus scaling argument x TS . The corresponding zero mode finite-size scaling 
functions (x, 0) are displayed as solid lines. Each data point was obtained by averaging over 10 5 MCS. 



magnetization histogram at the corresponding tempera- 
ture T C (L) and zero external field, i. e. at x rg = y rg = 0. 
Using these histograms, the considered quantities were 
extrapolated to nonzero fields an d plotted as implied 
by their finite-size scaling forms (see Fig. [31). The values 
of the fit parameters required for these plots are taken 
consistently as they were determined from the tempera- 
ture runs (cf. Table [i}. Since the reweighting technique 
allows the extrapolation of the quantities to arbitrary val- 
ues of the external field, the data are displayed as con- 
tinuous lines. 

Within the intervals of the scaling variables x rg and y rg 
that were considered in the simulations we find excellent 
agreement of the MC data with the zero mode theory. 
In Figs.EIa,b) merely for negative values of the tempera- 
ture scaling variable x rg remarkable deviations occur and 



the MC data with increasing system size L converge to 
the corresponding zero mode finite-size scaling function. 
This effect indicates the presence of further corrections 
that are expected to vanish in the limit L — > oo (see also 
Ref. that refers to the five-dimensional short-range 
nearest neighbor Ising model). 

The deviations visible in Fig. for large values of the 
field scaling variable y rg are due to the finite size of the 
histograms that were used to obtain the curves shown 
there by means of the reweighting technique. They orig- 
inate from the fact that a magnetization histogram of 
size K that is based on K values of the average spin s 
[Eq. JSJ] sampled in a MC simulation, in the mean does 
not contain values of s occurring in a simulation with 
a lower probability than 1/K. As will be demonstrated 
in the following, these deviations can also be reproduced 




Figure 3: Finite-size scaling plot of the magnetization m(t c (L) , h; L) (a) and m a b B (t c (L),h; L) (b), and the susceptibility 
x(t c {L), h; L) (c) and Xabs(tc(L), h; L) (d), versus scaling argument y rg . The size of the histograms that were used to obtain 
the curves amounts K = 2 x 10 6 MCS for each system length L. The zero mode finite-size scaling functions Y; mf (x, y; K) are 
plotted as solid lines where the cutoff parameter takes the value qa(2 x 10 6 ) ~ 3.2201 (see text). 



within the zero mode theory by truncating the zero mode 
probability distribution Vo(f), Eq. from which we 

start and perform the rescaling, Eq. I|24|). After that we 
set x m { — and y m { = since the magnetization his- 
tograms we used to obtain the curves shown in Fig. |3| 
were recorded in the simulations also for vanishing x Tg 
and y rg . The corresponding normalized zero mode prob- 
ability distribution denoted as Vo,c(<p) then reads 

VoAv) = (uL d )^l((uL d )i^) (40a) 
with the zero mode finite-size scaling function 

*^c(9) = J^-** 4 (40b) 



and hence decreases monotonically for increasing \ip\. So 
to reproduce the effect that is due to the finite size of the 
histograms we chop the tails of this distribution by lim- 
iting the integration range in Eq. (|17|l and the zero mode 
partition function, Eq. I]15p. to a finite interval [— <p\, (p\\. 
To make an appropriate definition of the cut-off parame- 
ter </?a = fiviK) > we use an estimation that is known 
from extreme value statistics [3(j • We assume that in the 
mean only one out of K measured values of the average 
order parameter tp lies outside the interval ] — tp\,tpA[, 
i.e. it fulfills \(p\ > tp a. This implies 

f°° 1 
P(M>^a)=2/ d^P ,c(v) = j7 (41) 
J VA K 

from which the cutoff ip\ is implicitly defined. Replac- 
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-11/2, -1/4 ., . 
% = C 2 SL ln (L/L q) 



Figure 4: Finite-size scaling plot of the critical histogram 
H c (s;L). The size of the histogram depicted here amounts 
K — 10 6 MCS for each system length L. According to 
Eq. this corresponds to a cutoff value of <2a(10 6 ) ~ 3.1728. 
The value of the fit parameter Lq and the nonuniversal metric 
factor C2 are taken as they are listed in Table Q 



the displayed histogram. The zero mode finite-size scal- 
ing functions Y^ ai {x, y; K) that were evaluated using the 
truncated zero mode probability distribution with the 
cutoff q\(K) are also shown in Fig. [3] and nicely agree 
with the extrapolated data. 

Before finishing this section we also want to discuss the 
importance of the corrections to the leading orders in the 
expansions in Eqs. Q that are determined by the param- 
eters Lq and v. For that purpose Figs. [S{a,b) show the 
finite-size scaling plot of the susceptibility Xabs{t,0; L) 
versus scaling variable x Tg (Fig.|^b)) where we set v = 
in Fig. |S{a) and Lq = 1 in Fig. GJb). The values of 
the other parameters entering the plots are taken from 
Table Q] As is expected, neglecting one of the cor- 
rections causes a significant displacement between the 
MC data and the zero mode finite-size scaling function 
y™^ a (o; r g, 0). This effect could partially be compensated 
by readjusting the remaining fit parameters, but this pro- 
cedure would in the case of Fig.[3Ja) lead to a wrong de- 
termination of the critical temperature T c (oo) (see also 
discussions in Ref. [3l| that refer to spin models above 
the upper critical dimensionality). 



V. CONCLUSIONS 



ing 'Po,c( ( /9) according to Eqs. 140fl the evaluation of the 
resulting integral yields 

r°° Ffl 27 4 \ -1 

2 £ dq Y?l (q) = ^ r( 2 f } gAj = 1, (42) 

where q\ = (uL d )i(pA and T(a, b) denotes the incomplete 
Gamma function. Expanding this expression and sepa- 
rating q\ = q\(K) yields to leading order a logarithmic 
dependence on the histogram size K of the form 

a -°° vf (i*®-^ 1 *®))* (43) 

with the constant K a = (3/4)*r(±). 

To demonstrate the validity of the definition of the cut- 
off </?a, respectively, qa, Fig. 0] shows a finite-size scaling 
plot of the critical magnetization histogram H c that was 
recorded for x Ig — y Ig — 0. Starting from the finite-size 
scaling form of the magnetization, Eq. (|9a|) . it is straight- 
forward to show that this histogram scales as 

H c (s;L) ~ C^I* ln-i(^)Y Hc (q Ig ) (44) 

with the finite-size scaling function Yh (q) and the scal- 
ing variable q rg = C^sL^ ln"(-^). As is expected also 
in Fig.0J the MC data are found to be in very good agree- 
ment with the zero mode expression (q) [Eq. (|40b() ] 
of the finite-size scaling function Yjj c (q) . It furthermore 
turns out that all data points lie accurately within the 
interval [— qj\(K), q^iK)] that corresponds to the size of 



We have introduced a two-dimensional long-range spin 
model that displays both isotropic and anisotropic phase 
transitions and, in particular, strongly anisotropic phase 
transitions. As a first stage the critical behavior of the 
model in the isotropic case for which it is found to be 
at its upper critical dimensionality was investigated. For 
that purpose we have carried out Monte Carlo simula- 
tions and studied the temperature and field dependence 
of several quantities. Using results of the renormalization 
group, the numerical data obtained for different system 
sizes were analyzed by means of a finite-size scaling anal- 
ysis. It turns out that beside a size-dependent shift that 
has already been discussed in the literature a character- 
istic length Lq that was inserted to the logarithms is an 
important correction that must not be neglected. 

Furthermore, the collapsed data were compared to the 
zero mode (mean-field) theory and found to be in excel- 
lent agreement. It turns out that the logarithmic cor- 
rections typically occurring at the upper critical dimen- 
sionality do only enter the finite-size scaling functions 
through their arguments, whereby these functions were 
derived from zero mode theory. This shows that at least 
in the present case the concept of universal finite-size 
scaling functions can be extended to the upper critical 
dimensionality. 

Finally we note that the numerical results strongly in- 
dicate the validity of the zero mode theory at the upper 
critical dimensionality and might shed new light on recent 
controversial discussions about its correctness for d > d u 

As it will be subject of a future work it is desirable 
to extend the analysis shown above to the anisotropic 
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Figure 5: Finite-size scaling plot of the susceptibility Xabs(i, 0; L) without the shift correction (v = 0) (a) and with Lo = 1 (b). 



case uiu 7^ wx- In particular, the critical behavior of 
the model should be investigated when approaching the 
strongly anisotropic cases lo\\ = —2uj±_ and — 2u/ii = c<j_l, 
respectively. 
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Appendix A: THE FUNCTIONS H m (x,0) AND 
H m (0,y) 

The evaluation of the function S m (x, 0) with Mathe- 
matica |33[ yields 

E m (x,0) = 3-i(™ +1 )4™(r(i2±i) lFl (^;l;x 2 ) 

- 2x r{^) ^(2^; §;**)) (Al) 

with the confluent hypergeometric function iFi(a;b;x). 
This expression is valid for x S K and m £ 1R >_1 and 
can be simplified further for a given integer value of the 
parameter m. For m = {0, 1, 2, 4} one obtains 

S o (x,0) = 3^e^Ti(i), (A2a) 

Si(x,0) = iy|e* 2 erfcOr), (A2b) 

E 2 (x,0) = -S-h^ftiW-xTiixj), (A2c) 

S 4 (x,0) = ^33 e #(2o: 2 T i (x)-3xT | (x)-T5(.T)) , 

(A2d) 



where the function 

T a (x) = ir(x 2 ) a (I~a(±x 2 ) - sgn(x)I a (±x 2 )) (A3) 

with the modified Bessel function of the first kind I a (x) 
has been introduced, and erfc(x) denotes the complemen- 
tary error function. The function T a (a;) is well behaved 
through zero argument for the pertinent noninteger pos- 
itive values of a. Assuming x > 0, Eq. (|A3|) simplifies 
to 



T a (x) x = 2x 2a K a (\x 2 ) sin(Tra) 



(A4) 



with the modified Bessel function of the second kind 
K a (x). 

An analogous treatment of the function ^ m (0,y) re- 
sults in 



Mo,y) = ~£(lf^&m (A .-„ 

k=l 



X 



r(fc) 

rp ( m+k i . k fc+1 fc+2 fc+3 . J/ 4 . 
2 4 \^ 4 ' ' 4 ' 4 ' 4 ' 4 ' 27 



where p F q (ai, . . . , a p ; b\, . . . , b g ; x) denotes the general- 
ized hypergeometric function. This expression is also 
valid for y S R and m S M >_1 and cannot be simpli- 
fied further for a given value of m involving less general 
functions. Finally let us note that the function E m (x, y) 
fulfills the recursion relations 



E m+2 (x,y) 



d 



2,^/idx 



^m(x,y), (A6a) 



d 



E m+ i(x,y) = —^ m (x,y) 



(A6b) 



as follows from Eq. I|26|) . 
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